Phase diagram of three-leg ladders at strong coupling along the rungs 
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A phase diagram of the t — J three-leg ladder as a function of hole dopping is derived in the 
limit where the coupling parameters along the rungs, t± and J±, are taken to be much larger than 
those along the legs, t|| and J|| At large exchange coupling along the rungs, J±/t± > -j-, there is a 
transition from a low-dopping Luttinger liquid phase into a Luther-Emery liquid at a critical hole 
concentration n cr u « 1/3. In the opposite case, J±/t± < -^=, there as a sequence of three Luttinger 
liquid phases (LLI, LLII and LLIII) as a function of hole dopping. 
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I. INTRODUCTION 

The recent experimental success in synthesizing quasi-one-dimensional ladder materials with a mobile charge carriers 
has raised an increased interest in the theoretical understanding of their rich phase diagram jij , || . In previous studies 
of microscopic models on various ladder geometries, a competition between superconducting, phase separation and 
density wave instabilities has been observed ||], (J]. In particular, it was seen that ladders with even and odd number 
of legs have quite distinct generic features, such as the presence of a spin gap at half-filling in even leg system 
and its absence for odd-leg ladders. When examining this half-filled case it was realized that inter-band scattering 
processes are relevant, and that this is the reason why approximations based on strong coupling anisotropies, such 
as expansions in J\\/J±, give a correct physical picture which extends beyond the isotropic regime. Furthermore, 
most ladder materials show coupling anisotropics within the ladder complex, e.g. a recent structural analysis of the 
vanadate ladder NaVa^Os suggests a strong rung-coupling anisotropy of J±/ J\\ ~ 4 ||. 
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FIG. 1. The t — J ladder with three legs and L rungs. The couplings along the legs are t\\ and J\\, and those along the 
rungs t± and J±. 



In this paper we further explore the strong rung-coupling limit of the t — J model in the presence of mobile holes 
. \ on a three- leg ladder (Fig.l). 

By comparing the results in this analytically tractable limit with numerical diagonalizations, we will see how far 
H \ this analysis can be extended towards the regime of isotropic coupling parameters. 

We arrive at a phase diagram valid in the limit J±,t± >> J\\, t\\. At Jj_ < ^t^J- h contains three various Luttinger 

Liquid (LL) phases of different nature, depending on the concentration of holes. For J± > -^t± at low hole doping, 
there is also a LL phase. However beyond a critical hole concentration a spin gap opens up, and a transition occurs 
to a Luttinger-Emery liquid (LEL) with an effective hole-hole attraction, as it has already been observed in previous 
studies of the isotropic case (Jj_ = J||, tj_ = t\\) ||, 0. 
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II. SINGLE-RUNG STATES 



The Hamiltonian of the anisotropic t — J model on a three-leg ladder is given by 

3 2 
H = "*ll E E P (ia^+1^ + H.C)P - t± J2 E P (4,a )CT <V+l, CT + B.c)P - 
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+J\\ ^ ^ fSi, a Si+i )0 — -n^ a n i+ i ia j + ( Jj_ ^ ^ Sj )0 Si ia +i — -n^ a n^ a+ i j 

i a=l ^ ' \ i a=l / 

where i runs over L rungs, c(=|J.) and a are spin and leg indices. The first two terms are the kinetic energy (P is 
projection operator which prohibits double occupancy) and the last two exchange couplings Jh (J±) act along the legs 
(rungs) . 

Let us start by discussing the low-energy states of H on a single 3-site rung with 0,1 and 2 holes in the limit where 
J|l = t\\ = 0. They are listed in Tab.l. In this limit, the exact ground state of the whole ladder is simply a product of 
these rung states. An important symmetry present in the 3-leg ladder is its reflection parity about the center leg, R. 
Along with the total spin quantum number, S, and its projection, S z , R characterized the symmetry of the ground 
state vector. 

Table 1. Ground state energies and vectors for 
a 3-site rung with 0, 1, 2, and 3 holes. 



E n 



4/1 



2tj_ 



V2tj 







ground state eigenvector 
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1000 > 



At half filling (0 holes) the ground state is two- fold degenerate (S=^; S z = ±|), and it has odd parity with respect 
to reflection about the center leg, R = — 1. As discussed previously [6,8], this state behaves as an effective spin 
-1/2 rung spin, and an inter-rung magnetic coupling Jn between such states introduces the low-energy behavior of a 
Heisennberg AFM spin -1/2 chain. 

Note that first excited state here is also (S — h; S z = ±^) doublet, but it has even reflection parity, R = 1, and 
thus is in a non-bonding configuration. A second excited state corresponds already to S = | and is irrelevant for our 
considerations. 

The ground state with one hole on a 3-site rung is a singlet (S = S z = 0) with even parity, R — 1. Subsequently, we 
will consider separately the two regimes of " strong coupling" Jj_ > > t± and " weak coupling" Jj_ < < t± . In strong 

coupling case E\ w — Jj_ — 2^, while in weak coupling case E\ i=s — 4^ — \/2t_\_. It will be shown later that hole-hole 
pairing on a rung (leading to LEL) occurs naturally in the strong coupling limit, while it is absent at weak couplings. 

The first excited state corresponds to a nonbonding singlet with R = —1, while a second excited state - to an 
antibonding singlet with R= 1. 



For 2 holes on a 3-site rung, the ground state is a (S 



S z = ±i) doublet with R = 1. The first excited state 



is non-bonding with R = — 1, while a second excited state is antibonding with R = 1 again. 

So let us compare now in the limit of almost independent rungs (Jn = tu = 0) the ground state energies for different 
configurations. As a result we obtain: 



1. for hole concentration < n < i a minimal energy minE = E a = (1 — 3n)EoL 



inE\L -corresponds to a 



mixture of rungs with one hole and without holes (Fig. 2) 



FIG. 2. Configuration a). Rungs with one hole in the surrounding of rungs without holes 
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2. for | < n < | - two configurations are possible: 



E b = (2- 3n)EiL + 3nE 2 L 
E c = |(1 - n)E x L + (|n - 



(1 - n)£i£ 



(2) 
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FIG. 3. Configuration b). Rungs with one hole in the surrounding of rungs with two holes. 
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FIG. 4. Configuration c). Rungs with one hole in the surrounding of rungs with three holes. 

Configuration b) is given by Fig. 2 and corresponds to the mixture of rungs with one and two holes, while 
configuration c) is given by Fig. 4 and corresponds to the mixture of one and three holes. 

3. Finally for | < n < 1 there are two possibilities again: E c = |(1 — n)E\L - is given by Fig. 4 again while 
Ed = 3(1 — n)E2L — corresponds to a mixture of rungs with 3 and 2 holes (Fig. 5 ) 



FIG. 5. Configuration d). Rungs with two holes in the surrounding of rungs with three holes. 

Comparison of Ed and E c yields: E<x = \E\ and hence: —2^/2t± = 
As a result: 

i^L) - JL 

[ t± )crlt ~ V2 



(3) 



For J|_ > J± C rit in ah the region | < n < 1 configuration c) [3+1] holes is realized. ForJj_ < J± C rit in the 
region | < n < | configuration b) [1+2] holes is realized, while for | < n < 1 configuration d) [3+2] holes is 
more beneficial, (see Fig. 6) 
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FIG. 6. Phase diagram for different ground state configurations as a function of J±/t_i 



III. KINETIC ENERGY OF RUNGS DELOCALIZATION 

Let us " switch on" the next approximation and calculate kinetic energy of rungs derealization. 
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To be more specific for [0+1] phase we need to calculate a matrix element which interchange the rung with holes 
and the rung with 1 hole. (Fig. 7) 



1. 

a 2< 
3- 

i - 1 / i - 1 i 

FIG. 7. Interchange of the rung with holes and 1 hole in first order in tii. 

On electronic language we calculate: 

~t\\ E c ta^-ha,a |*<(2e0*i-i(3el)) = -t|| ^ |< a)CT *i(2eZ)> Ic^^^^eO) , 

a, cr a,cr 

where c^~ a CT , Ci_i i0 , CT are electronic operators; 

^i(2el) = ^i(lh); *i_i(3eZ) = *i_i(0/i) (see Table 1) 

Note that large scales J_l and i_L fix the global configuration, i.e. for given concentration of holes both the number 
of rungs with 2 electrons and the number of rungs with 3 electrons do not change. 
That is why: 

-t\\ |<a, CT **( 2e 0) \ci-i, a ,a^i-i^el)) = -t e// *i(3el)*i_i(2el). 

a,a 

The calculation of the i e // yields: 



3i|| a 2 

4 Q!2 — Ot\ ' 



(4) 



where cti^ = — 2tt ~ ~ C0Cmc i crL ts which enter, respectively, in the eigenvectors of a ground state (see eq.(3)) 

and an antibonding state for 2 electrons on the rung. 

Note that to calculate t e ff we used anticommutation relations for fermionic operators c + and c together with a 
condition which prohibits a double occupancy. 

In the limiting cases expression (4) reads: 



l eff 



In the language of effective operators: 



3t|j 
4 

3*1 



1 + 

1 



3/l 



iJ±arit 



for Jj_ >> J± C rit, 
for Jj_ << J± C rit- 

T a i-l,a- 



(5) 



An operator a^" CT corresponds to the creation on the site z of a rung with 3 electrons and simultaneous destruction 
on the same site of a rung with 2 electrons. So, we could represent a^ a in the following form 



(6) 



Here an operator f?~ a creates the rung with 3 electrons and total spin S — 1/2. Hence it has a fermionic nature. 
An operator bi destroys 2 electrons (an electronic singlet with total spin S = 0) and hence has a bosonic nature. 
Of course, 2 rungs could not occupy the same place. It means that they are subject of infinitely strong Hubbard 
repulsion: 



Uoo ^ ^ Pi, a Pi, — a > 
la 



(7) 
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where p ia = af a a iia = (f^ a fi tlJ )(b~l a b itCr ). So, in a state [0+1] holes a system is described by effective Hamiltonian: 

Ho+i = -Efl+l — t e ff a t,a a i,<? + Pi,aPi,-a- (8) 

icr i 

It is ID fcrmionic Hubbard model with repulsion. We know that it belongs to the universality class of Luttinger 
liquid §. 

/ \ 1/2 

Note that a more detailed analysis in case of J± >> t± shows that at densities n* > ( J (which, in principle, 

could be smaller than 1/3) we will have a two-band degenerate Hubbard model instead of a one band model. However, 
this situation will also fall in the universality class of Luttinger liquid I. 
Let us consider now a state [2+1]. Here: 

t eff = (* i (2el)* < -i(le0 Hiic^Ci-i.^l ^(leZ)^_i(2eZ)) 

Direct calculation of i e // yields: 

tu 1 



2 (a 2 - ai) 



For the limiting cases: 



1 + 2i 



5/ 



3J ± 
3Ji 



20J, 



«i \ (\/2 + ai) 
2^ J ° 2 2~ 



for Jj^ >> J± cr it, 

for Jj_ << Ji crit . 



(9) 



As a result: 



An effective operator = !ta°i nas a fermionic nature again (see Fig. 8). It creates a rung with 1 electron and 
5 = 1/2 and simultaneously destroys a rung with 2 electrons and 5 = (bi). Finally: 



where 



#2+1 — E 2 +l — t e ff 2^ d^ a di-\ t a + C^oo Pi,aPi-a, 



Pi, a — a di,a 



(10) 
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FIG. 8. Interchange of the rung with 1 electron and the rung with 2 electrons due to tu. 

We again derive a ID fermionic Hubbard model with repulsion. So a state [2+1] corresponds to LLII which describes 
a motion of a rung with 1 electron in the surrounding of rungs with 2 electrons. 
Now let us proceed to the case [3+2] (Fig. 9). Here 

teff = (%(lel) |-t|| C £ Cj_i )0)ff | *i_i(leZ)> 




/ - 1 i i-\ i 

FIG. 9. Interchange of the rung with 1 electron and an empty rung due to tu. 
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It is easy to derive that t e // = t because this problem is equivalent to the motion of an electron in the empty space. 
Hence: 



Hkin - -t e ff ^2 



i—l.ai 



where fl~ = ^(leZ) - has a fermionic nature again. (It creates a rung with 1 electron and destroys an empty rung). 
As a result: 



H- 



3+2 



(11) 



where Xi >a = f^ a fi irr is a density of rungs with one electron. 

Hamiltonian (11) again corresponds to ID Hubbard model with repulsion. So, it describes LLIII where the rungs 
with 1 electron move in the surrounding of rungs without electrons. 

A state [3+1] has quite a different nature. 

t 2 

Here t eff ~ r^=^\ (see Fig. 10). 



t, 



FIG. 10. Interchange of the rung with 2 electrons and an empty rung in second order of perturbation theory in t||. 



Since a phase [3+1] is realized only for J± > J 



J-crit- teff 



of the rung with 2 electrons in the surrounding of empty rungs. Hence 

Hkin = -teff ^ b t b i-i> 



ijj/ J±. Kinetic energy in this state describes a motion 



(12) 



where bf = 1 $i(2el) is a bosonic operator which creates a rung with 2 electrons and destroys an empty rung. Of 
course, 2 rungs with 2 electrons can not occupy the same place. As a result: 



#3+1 — E3+I — teff ^ bfbi-1 + Uoc ^ 



n bos n bosi 



(13) 



= bfbi, and now we have ID Bose- Hubbard with strong repulsion. This model belongs to universality class of 
Luther-Emery liquid. It has a spin gap at half filling and large superconductive fluctuations in a doped case. When 
we include a boson rescattering between neighbouring ladders, than a finite Tq arises in the system p0| . 
Finally, the phase diagram reads (Fig.ll): 
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LLII 
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FIG. 11. The phase diagram of three-leg ladder at strong coupling along the rungs. 

IV. THE ROLE OF AFM EXCHANGE ALONG THE LEGS J, . 

In a close analogy with a double exchange model for JhS >> zt > JffS 2 : 



G 



- a largest scale (FM exchange Jh) forms a local onsite state with Stot = S + 1/2 and then effectively drops out of 
the model. Low energy physics (including phase separation on FM and AFM regions) is governed solely by smaller 
parameters t and Jf f . Absolutely the same scenario is realized in our model. The largest parameters J± and t form 
the stable configurations LLI, LLII. LLIII, and LEL and after that effectively drop out of the model. Low-energy 
physics is governed solely by Jm and in. 

Then, by analogy with FM-polarons formation in double exchange model, we could have in our case either a diluted 
configuration (Fig. 12) or a phase-separated state (clusterization) (Fig. 13). The clusterized phase was found in in 
numerical study of an isotropic regime i|| = J|| = Jj_; J = 0.35t. 




FIG. 12. A diluted configuration in LLI 



FIG. 13. Clusterization in LLI 

A diluted phase is more beneficial with respect to kinetic energy ~ while a clusterized phase is more beneficial 
with respect to magnetic exchange energy between the rungs ~ Jm. So, in the case of small J|| < t\\ a diluted phase 
is more beneficial. The energy of this state in case of LLI reads: 



E Q+1 = (1 - 3n){E + Eq Q )L + 3n(E 1 + E^ n - E J m )L. 



(14) 



Here Eq™ = —2t e fj (where t e fj is given by (9)) - corresponds to derealization energy when rungs with one hole 



occupy the bottom of the band 

< - (*, + (3eO*,--i(3eQ 



J\\ ^ ^ (S^qSj^l^ — — ^i^^l-l.o) 



*+(3eZ)*+i(3ei) 



-corres 



;ponds to a chessboard AFM configuration of rungs with holes (see Fig. 14). 



© O © 
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FIG. 14. A chessboard AFM configuration of rungs with holes. 
Effective Hamiltonian for LLI reads: 

H = Eq + \ — t e ff ^ af a ai-i tCr + Uoo ^ Pi,aPi-\,a + J\\ ^ SiSi_i, 



(15) 



where Si is a rung spin on site i. 

So, a final effective H is given by i e // — J\\ model. 

For J|| < 2t e ff it still belongs to the universality class of LL. There is no spin gap at half- filling. The basic instability 
for moderate doping is towards SDW- formation |ll| . 

For J|| > 2t e ff, in total analogy with 2D t — J model, a bound state appears and model becomes unstable towards 
superconductivity. So, for small J|| < i e // its role is just: 

1. to form small energy corrections, 
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2. to change a little bit phase boundaries on a phase diagram: 



1 



1 



1/2 



?y is a numerical coefficient, 
3. to introduce in the model effective AFM attraction between neighbouring rungs. 
As a result the gross features of phase diagram on Fig. 11 are conserved. 

V. DISCUSSION 



Now let us proceed from the calculations to qualitative arguments. In isotropic case, besides a tendency towards 
clusterization |7|, there is a tendency towards coexistence of LL and LEL. (Fermi-Bose liquid scenario restated recently 
in (TJ]). In strong coupling regime there is a gap between the bottoms of Fermi-gas band and a Bose-gas band. The 
gap is of the order of t± — 2t\\. In an isotropic regime we could overcome this gap due to an increase of tii (in = tj_) 
and obtain a picture of Fig. 15: 



W22E!m~~~ TjmmmA 

Bose-gas 

Ferm,-gas 

(LL) ( ; 

FIG. 15. Fermi-Bose liquid scenario restated recently in [12]. 



VI. CONCLUSION 



In conclusion let us emphasize that 

1. 3-leg ladder possesses both properties of 2-leg ladders and ID doped spin chains. 

2. There is a hope that when we increase a number of legs the difference between odd and even numbers will 
become smaller. In favour of this assumption is a fact that in 2n-leg ladders spin-gap scales as: 

A 2 „ ~ A 2 /(2") - Jl/2™ for n -> oo. 



3. HTSC-materials have both properties of two and three leg ladders. 

At n < (2 -j- 4)% (AFM region) there is no spin gap in HTSC materials as in 3-leg ladders for J± <t±. 

However at (2 4)% < n < 20% there is a spin pseudogap in HTSC materials in analogy with 2-leg ladders. It 
means that HTSC has a long prehistory when we go from half-filling (n = 0) to an optimal doping values (n — 0.15). 
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